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Аппроксимационное построение математических моделей по точечным 
экспериментальным данным методом си-дше” 


Р. А. Нейдорф 


Рассматривается решение проблемы описания экспериментально полученных зависимостей. Предложенный 
автором метод основывается на получении локальных приближений фрагментов данных зависимостей и ад- 
дитивном сведении их в единое аналитическое выражение. Такой эффект определяется применением спе- 
циальных «выделяющих» функций, ограничивающих область ненулевого определения каждой из аппрокси- 
мирующих функций. По применяемым принципам аппроксимации метод назван «сиЁд/ие»’. Наиболее близ- 
кий его аналог — сплайн-аппроксимация. Однако метод сиё -д/ие значительно гибче, т. к. не связан ни коли- 
чеством аппроксимируемых сплайном точек, ни порядком функций, аппроксимирующих участки. Порядок 
аппроксимирующего полинома или другой аппроксимирующей функции, а также её структура для каждого 
участка могут быть произвольными. Ещё одно выгодное отличие сиЁо/ие аппроксимации состоит в единой 
аналитической записи всей кусочной функции вместо задания векторной сплайн-функции громоздкой систе- 
мой уравнений. Такого эффекта удалось добиться применением аналитической функции, аппроксимирующей 
и параметрически сколь угодно приближающейся к единичной функции Хевисайда. Приведены результаты 
аналитического и численного исследования свойств и итогов применения предложенного метода. Они про- 
иллюстрированы практическими примерами приложения метода к практическим задачам, табличными и 
графическими данными. 
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тивность, аддитивность, дифференцируемость, аналитическая функция, параметрическое приближение. 


Введение. Одним из важных практических приложений математических методов является реше- 
ние задачи аппроксимации различного рода экспериментальных зависимостей. К аппроксимации 
приходится прибегать при построении математических моделей различных природных, техниче- 
ских, экономических, социальных и иных явлений и объектов на основе результатов эксперимен- 
тальных исследований [1, 2]. Наиболее распространённым и математически проработанным под- 
ходом к решению таких задач является регрессионный анализ, опирающийся на метод наимень- 
ших квадратов (МНК) [1, 3]. 

Метод регрессионного полиномиального описания экспериментальных точечных данных в 
заданном диапазоне исследования является стандартным. Он поддерживается всеми универсаль- 
ными и несколькими специализированными программными средствами (МАТЁЬАВ [4, 5], МайСАО 
[6, 7], ЗайНса [8, 9] и др.). Возможны и другие, не полиномиальные формы регрессионной ап- 
проксимации, когда эмпирически подбирается совокупность подходящих для описания функций 
[7, 9]. Но этот подход нерегулярный и не поддержан таким мощным, хорошо разработанным ма- 
тематическим аппаратом, как классический регрессионный анализ. Однако и у последнего есть 
существенный недостаток — он плохо подходит для аппроксимации многоэкстремальных зависи- 
мостей, особенно кусочного характера с явными изломами, т. е. с фактически существующими 
разрывами по производной. В этом случае применение при обработке данных критерия МНК мо- 
жет дать два различных результата. Первый, при малом порядке аппроксимирующего полино- 
ма — очень грубое приближение к экспериментально снятым точкам. Второй, при высоком поряд- 
ке — достаточно близкое воспроизведение зависимости в окрестности этих точек. При этом не 
обеспечивается гладкость описания межточечных интервалов. Так, на рис. 1 приведены примеры 


* Работа выполнена по темплану Минобрнауки России № 2.3.13 «Метод сбалансированной дискретизации для задач ими- 
тационного моделирования динамических процессов в распределённых объектов». 
* Резать-клеить (англ.). 
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результатов регрессионного описания экспериментальной точечной зависимости от угла наклона 
а аэродинамического коэффициента подъёмной силы (рис. 1) одной из конструкций аэростатов 
[10-12]. Сила возникает либо при воздействии бокового ветрового потока (а — угол крена), либо 
при движении аппарата (а — угол атаки). Зависимость необходимо оценить по 13 эксперимен- 
тально полученным точкам. В связи с этим максимально возможная степень оптимального по МНК 
аппроксимирующего полинома равна 12. Регрессия этого порядка построена на рис. 1, а. Хорошо 
видно, что график полинома 12 порядка абсолютно точно воспроизводит экспериментальные точ- 
ки, но практически произвольно ведёт себя на межточечных интервалах. 
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Рис. 1. Примеры аппроксимации точечной экспериментальной зависимости от порядка аппроксимирующего её полинома, 
оптимального по МНК: а— п= 12; 6— п=9; вл п= 6 


Это полностью согласуется с известным свойством результатов классического регрессион- 
ного описания — с повышением порядка оптимальный по МНК полином теряет фильтрующие 
свойства [1]. На нижних графиках рис. 1, 6, в показаны результаты регрессионного описания ис- 
следуемого диапазона полиномами девятого и шестого порядков. Видно, что аппроксимация 
межточечных интервалов последовательно улучшается, но воспроизведение экспериментальных 
значений аэродинамического коэффициента (АДК) резко ухудшается. Таким образом, полноин- 
тервальная полиномиальная аппроксимация моделируемой зависимости, имеющей очевидный 
кусочный характер, в общем случае неэффективна. 
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Для подобных зависимостей более продуктивен другой подход, основанный на математи- 
ческом описании отдельных интервалов кривой и записи зависимости в виде системы функций со 
смежными областями определения каждой из них (кусочная аппроксимация [13, 14]). Такая 
функция записывается в виде логически обусловленной по областям определения системы функ- 
ций. Это делает результат описания громоздким и не очень удобным для компьютерного модели- 
рования. Кроме того, при использовании функции дробной структуры дополнительные проверки 
логических условий по значению аргумента и обращение к новым расчётным формулам создают 
тормозящий эффект при онлайн-моделировании быстро протекающих процессов. Дополнитель- 
ные трудности, вызванные кусочностью, возникают также при необходимости использования 
производных от моделируемой зависимости. 

Более совершенным является метод сплайнирования, или сплайн-функций. Его отличи- 
тельная особенность — совмещение кусочной аппроксимации с краевым согласованием значений 
производных. Последнее делает функцию дифференцируемой, несмотря на кусочность [14, 15]. 
Особенно распространены кубическая сплайн-аппроксимация и интерполяция, обладающие высо- 
кой регулярностью решения и хорошими точностными характеристиками [16]. Однако высокая 
степень дискретизации, жёсткая привязка к количеству сплайнируемых точек и громоздкость си- 
стемы условно логического описания делают этот подход неудобным для использования при 
необходимости построен 

ия компактной математической модели для аналитических преобразований и имитацион- 
ного моделирования. Кроме того, основное достоинство сплайн-подхода — согласование аппрок- 
симирующих функций по производным не всегда физически оправдано для некоторых объектов. 

Таким образом, очевидна актуальность поиска и исследования новых подходов, методов и 

алгоритмов экспериментально-аналитического описания зависимостей, имеющих очевидный ку- 
сочный характер, указывающих на явное наличие разрывов первого рода по производной на гра- 
ницах участков. 
Постановка задачи. Ставится задача построения метода аппроксимации, позволяющего созда- 
вать математические модели зависимостей, представленных экспериментально полученными то- 
чечными данными, в виде единой аналитической, хотя бы один раз дифференцируемой функции. 
Перечисленные свойства позволят строить аппроксимационные математические модели любых 
зависимостей, включая кусочные. 


Общий подход к решению задачи. Пусть исследуемая зависимость РЁ (ух) задана множе- 


ством // точек (хе, у) в пространстве А?, где 1 =1, №, а верхний индекс е отмечает экспери- 


ментальный способ получения величины. Пусть это множество может быть разбито на л непере- 
секающихся подмножеств 2. (у,х). То есть У/,К=[1,п |-> (у,х) ПА, (у,х)=@. При этом 


1 


С. 


И 
СА 


должно выполняться условие |_/Е. (у, х)=Е (у, х). 


1 
Рассмотрим основополагающую идею решения поставленной задачи. В качестве первой её 
составляющей сформулирован подход, основывающийся на возможности аналитического описа- 


ния с заданной точностью 6%? любого выделенного /-го участка аппроксимируемой точечной 


зависимости в выбранной локальной области определения х е[Ху1х, | некоторой функцией 


Ф, а Данная функция обладает следующими свойствами: 


с 


У; = ф, (х,): Ух; [хиих, > УЕ -У; 





«5%, (1) 


где 6%? — допустимая погрешность аппроксимации на /-м участке; верхний индекс с означает 


вычислительный способ получения величины. 
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Метод и качество аппроксимации представляют собой отдельную задачу и в данном ис- 
следовании не рассматриваются. В общем случае эта задача решается для каждого участка мини- 
мизацией некоторого функционала О, [Е (у,х), Ф, (х)], формализующего оценку точности ап- 
проксимации каждого участка [14]. Подобный подход значительно повышает точность аппрокси- 
мации каждого участка и, как следствие, общую точность аппроксимации всей моделируемой за- 
висимости. Один из возможных и эффективных теоретических инструментов реализации такого 
подхода — регрессионный анализ, а наиболее распространённым и эффективным критерием ми- 
нимизации ошибки моделирования О, | * | является квадратичный критерий". 


Вторая составляющая предлагаемого метода решения — предположение, что можно скон- 
струировать некоторую функцию # (х) так, чтобы во всем остальном диапазоне оси х, кроме 


описываемого участка [ых х, | ‚ она имела практически нулевое значение. То есть 


УХ = [Х,4их, |-> (х)=Ф, (х). (2, а) 
В остальной области определения должно выполняться неравенство 
Ух: х,.- Ах, >х>х, + Ах, >Ё(Х)=0. (2, 6) 


Например, можно задаться условием 
Ухе[х,, - Ах, „Хх, +Ах, | (х) < 5, 
где Дх,, 077 >0 выбираются достаточно малыми в соответствии с точностными требованиями 
решаемой задачи. 

Функция со свойствами 2, а, 6 названа в этой работе «интервально выделенной» (ИВФ). 
Она как бы реализует операцию сиЁ — вырезания из всей точечной зависимости хорошо аппрок- 
симируемого участка. 

Располагая ИВФ для каждого обособленного по конфигурации интервала аппроксимируе- 
мой точечной зависимости, можно построить обобщённую аппроксимирующую функцию 


п 
РЦ.) (3) 
та 
Она будет описывать моделируемую кусочную зависимость с точностью 
б„ = тах {5,| /=[1,п]}, 
где б, =0%7 +6>?, 6, — совокупная ошибка аппроксимации функцией (1) каждого / интервала, 
полученная при построении кусочной функции Ё (х); б, — наибольшая ошибка, возникаю- 


щая при построении кусочной функции / (х). 


Выражение (3) представляет операцию д/ие — как бы склеивания всей точечной зависи- 
мости из совокупности ИВФ. В результате Ё (х) (3) будет единой аналитической функцией неза- 


висимо от количества порознь аппроксимируемых интервалов и их размеров. Эта функция ап- 
проксимирует весь описываемый точечной (в том числе кусочной) зависимостью диапазон и реа- 
лизует таким образом заявленный в постановке задачи метод сиЕод/ие. 

Как явствует из приведённых выше положений, сиЕо/ие аппроксимация реализуется в три 
этапа. На первом экспериментально полученная точечная зависимость разбивается на обособ- 
ленные интервалы и осуществляется математическое описание каждого участка какими-либо 
аналитическими функциями ф, (х) со свойствами (1)", хорошо их аппроксимирующими. На вто- 


* Соответствующий метод решения регрессионных задач известен как метод наименьших квадратов (МНК). 
* Этот этап удовлетворительно решается известными методами [1, 3, 13, 14, 17], в частности, регрессионного анализа, но 
такая задача в данной работе не рассматривается. 
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ром этапе для каждого участка формируется интервально выделенная функция / (х) ‚ обладаю- 
щая свойствами (2). Это позволяет на третьем этапе получить результирующую обобщённую ап- 
проксимирующую функцию (3) в аддитивной форме. 

Мультипликативная реализация интервально выделенной функции. Суть предлагаемого 
метода состоит в том, что для формирования ИВФ используется специально сконструированная 


«мультипликативно выделяющая интервал» функция (МВФ) — ^(х, 5 } Её задача — сужение 
фактической области определения выбранной для аппроксимации точечного подмножества / -го 
интервала аналитической функции ф, (х), которой обычно является вся числовая ось. Такое 
сужение достигается путём перемножения ИВФ и МВФ: 

Е (х)=Ф, (х)-Л(х, 5, (4) 
что и определило использование термина «мультипликативная». 

Сужение необходимо произвести до подобласти, соответствующей моделируемому интер- 
валу [х,.‚‚х, |. Поэтому в диапазоне значений аргумента [х,.‚х,| функция ^(х, 5,) должна 
принимать единичное значение. При всех других значениях аргумента МВФ (ввиду мультиплика- 
тивности её применения) должна принимать значения, сколь угодно близкие к нулю”. Тогда 
К (х) в (4) принципиально может отвечать свойствам (2, а, 6). Для достижения этого результата 
в МВФ введён второй аргумент $, функции ^(х, 5; Он подразумевает множество математиче- 
ских «настроечных» параметров $, = 5,, которые обеспечивают функции, необходимые для ре- 
шения задачи свойства. Имеется в виду, что варьирование 5. В МВФ должно обеспечивать до- 


стижение свойств (2, а, 6). 
Для решения поставленной задачи в данной статье предлагается вариант МВФ, обеспечи- 
вающей указанные свойства и делающей возможной сиЕо/ие аппроксимацию. Решение задачи 


рассматривается применительно к обработке одной функции Ф(х), аппроксимирующей точки 
произвольно выбранного интервала | х,,х, | ', но не совпадающей со всей кусочной зависимостью. 


Пусть искомая функция Л(х, 5,) задана следующим выражением‘: 
|х-х, + (х-х,) +2 х, -х+\/(х, 5% += 


4. [< -х,) + = [(х, -х) + =” 
В функцию (5) включено минимально необходимое множество параметрических настроек МВИ- 
функции: 








(5) 


Меж) = 


ЕЕ: 
Здесь х,, х, — левый и правый края аппроксимируемого диапазона; = — показатель погрешно- 
сти аппроксимации. Координатные настойки х, и х, обеспечивают формирование области 
[х„›х, ‚а = — выполнение условий (2, а, 6). 


Исследование функции (5) как многопараметрической зависимости показывает, что она 
обладает шестью чрезвычайно важными для решаемой задачи свойствами. 


* Степень близости определяет точность аппроксимации. 

* Применительно к рассматриваемому произвольно ориентированному на х интервалу / =/-1 — левый край его диапа- 
зона, г =/ — правый. 

* Функция построена на основе принципов, разработанных автором для решения задач квазиоптимизации быстродействия 
законов управления [18] и обнаруживших при исследовании более широкий диапазон возможных применений [2]. 
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1. Функция (5) симметрична относительно средины диапазона | х,,х, |. 


Действительно, подстановка в (5) симметричных относительно центра диапазона 











(х.+х,) (х,+х,) (х.+х,) . 
х=—^>_ аргументов, например х =" + бих 707 даёт одно и то же вы- 
ражение 
2 2 и 2 з 
х.=х х.-х х.-х Х=Х 
0,25. Жо 5 +=? |. в) 9,5 +=| . (< р Ее) и 
2 2 2 2 
зависящее от констант функции ^ (5) — х,,х,, = и смещения б. 
, (х.+х,) В 
2. Функция (5) имеет единственный максимум в точке хт =-—^ ——^ , соответствующей 
средине диапазона [х,,х, | 
0,25 
Х(хт, х,Х,, =) =тахл(х,...) = Лт = 0,25 + я (6) 
1+ е- 5 
(х„-х,) 





и два нулевых инфимума, соответствующих пределам х > —ю и х-›®. 
Иными словами 
№ Л(х,Х,,Х,, =) =0; 


х>-— 7 
т Л (х, хи х,, =) = 0. (7 


Хо 


Действительно, вычисление частной производной по х от (5) для нахождения необходимых 
условий экстремума даёт выражение 


2% (5) о - [© -х} +=] Не 











3 


ОХ 4-ю +=], -х}+е] (8) 


+3 (ж? +х,2)+2- =? |+ х-[ж+х? + (х, +х,)-=? ] 





(х.+х,) 


Полином в числителе этой функции обращается В НОЛЬ В ТОЧКе ХТ = и асимптотически 


стремится к нулю при х > +. 
В точке хи? знаменатель дроби (8) не нулевой. Его значение определяется полиномом от 
2 
настроечных параметров функции (5) (4-2 - (х, -х,) }> О и при любых х,,х,, = положитель- 
но. Значения знаменателя (8) при асимптотическом стремлении х -> + стремятся к +ю. 
3. Функция (5) в диапазоне [2 х, | может принимать значение, сколь угодно близкое к 


единице. Причем как близость к единице, так и степень охвата диапазона эффективно настраи- 
ваются параметром =. 
Действительно, из анализа выражения (5) следует вполне ожидаемый результат, показы- 
вающий что 
ит (ЖИТ, Хи, Х,, =) =1 


при любых х,,х,. 
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Кроме того, можно показать, что 
УХ, Хи ХХ, < Хи ХЕ [хи х,[ > ТА (Х, Хи Х,, Е) =1. (9) 


Таким образом, значения функции ^(х, Жк Е) могут быть приближены к единице с любой 


требуемой для задачи точностью во всех точках интервала, исключая краевые. 

4. Функция (5) в диапазонах [ -о, х, | и | х,, ®] может принимать значение, сколь угодно 
близкое к нулю. Причём и близость к нулю, и близость к краевым точкам значения функции так- 
же эффективно настраиваются параметром =. 

Действительно, анализ выражения (5) показывает, что эта функция асимптотически при- 
ближается к нулю по мере удаления аргумента от выделяемого диапазона, т. е. 

№ (ХХХ, ,=)=0 


Хх—>—5 
Хх 
при любых х,,х,,Е. 
Кроме того, можно показать, что 


УХ, ХИ ХЕХ, < их [Хих, | > ТА (Х, Хи Хх, Е) =0. (10) 
Таким образом, значения функции ^(х, Хи Х,, =) могут быть приближены к нулю с любой требу- 


емой для задачи точностью во всех точках вне границ интервала, исключая краевые. 
5. Функция (5) в краевых точках принимает значения 
Ух хх, <х, > А(Ххи Хи Хх, Е) =А(Хи ХХ, Е) = 0,5 (11) 
при любых =. 

Это легко проверяется прямой подстановкой. 

6. Функция (5) бесконечное число раз непрерывно дифференцируема, как и любая функ- 
ция с дробным показателем степени. Это свойство чрезвычайно важно для использования ап- 
проксимационной математической модели в задачах аналитического исследования. 
Графическая иллюстрация свойств функции мультипликативного выделения интер- 
вала. На рис. 2 приведён пример, иллюстрирующий вид МВИ-функции для настроечных парамет- 
ров х, =-4, х, =4, &==0,3. Величина = завышена, чтобы показать аналитичность поведения 





функции вблизи границ выделяемого интервала. 


2 





^(х,хььх.:0,3), „(х,хььх:»0,3), ^.(х,х1»Х:›0.01) 
































Рис. 2. Графическая иллюстрация поведения МВФ при различных = 
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Здесь же показан вид этой функции, приближённый к прямоугольному импульсу при 
= = 0,01. При меньших = плавные переходы на изломе графика не различимы. Таким образом, 
параметр = является эффективным инструментом настройки точности аппроксимации. 

На рис. 2 представлен также график производной МВФ. При уменьшении = импульсы по 
краям интервала приближаются к дб -функциям. 

Выявленные шесть свойств МВФ показывают возможность реализации интервально выде- 
ленной функции на основе мультипликативной модели вида (4), которая качественно отвечает 
требованиям к функции Ё (х) ‚ сформулированным в виде условий (1—3) в общем подходе к ре- 


1 
шению задачи. 
Однако наряду с качественной возможностью реализовать «вырезание» участка функции 


в виде Ё(х) из полной зависимости ф,(х) необходимо обеспечить количественную сторону 
условий (1—3), т. е. необходимую близость этих функций внутри интервала [%Х, | и близость 


Е (х) к нулю за его границами. Аналитическое исследование и решение такой задачи является 


р 
предметом достаточно серьёзных и громоздких математических преобразований. Оно выходит за 
рамки данной работы. 

Поэтому в статье приводится частная количественная оценка влияния = на точность ин- 
тервального «вырезания», реализуемого функцией (5). Для иллюстрации этого эффекта на рис. 3 
построены графики зависимости значения функции Л(х, -Л, х,,х,,=) от = для ряда различных 


А =0,005; 0,01;0,02; 0,04; 0,1. 


0.5 





Т 
^.(Х-А хх») 








0.4 








0.3 
































0.1 





















ра > 
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Рис. 3. Иллюстрация зависимости погрешности «вырезания» МВИ-функцией заданного интервала [х,,х, | от = при раз- 


личных Д. Сверху вниз кривым соответствуют значения Л: 0,005; 0,01; 0,02; 0,04; 0,1 


Видно, что с увеличением =, откладываемого по оси абсцисс, и с уменьшением Д, опре- 
деляющего близость исследуемой точки к границе интервала (кривые, располагающиеся выше по 
оси ординат), точность реализации нулевого значения уменьшается. Асимптотическим пределом 
этой ошибки является значение МВФ на границе, т. е. 0,5. 


Характеристики точности реализации нулевого значения Л(х, +Д,х,,х,,=) справа от 


правой границы х, зеркально аналогичны приведённым на рис. 3 в силу свойства 1 симметрич- 
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ности функции (5). Такую же оценку точности реализации единичного значения можно получить 
для точек ^(х,+А, х, Хх, Е) и Л(х, -А, Хх, х,, Е), находящихся внутри интервала, т. к. относи- 


тельно точек пересечения абсцисс х, и х, левая и правая ветви МВИ-функции при бесконечно 


большом интервале симметричны. На конечном интервале встречные ветви составляющих внут- 
риинтервального значения функции несколько снижают точность реализации единичного значе- 
ния при конкретном =, но ошибка может быть сколь угодно уменьшена снижением величины этого 
параметра. 

В поле основного графика дан наиболее полезный для оценки и выбора нужной величины 
= фрагмент начального поведения кривых. При этом параметр Д выбирают из допустимой вели- 
чины «размытия» выделяемого интервала, а = — из допустимой точности реализации как единич- 
ного значения внутри интервала, так и нулевого значения за его пределами. 

Пример реализации интервально выделенных функций и сиё-д/ие аппроксимации ими 
кусочной точечной зависимости. В качестве примера эффективного применения разработан- 
ного метода рассматривается задача математического описания зависимости одного из парамет- 
ров т. н. «стандартной атмосферы»” — температуры 7 — от высоты И. Точечные данные из таб- 
лицы, приведённой в [19], построены на рис. 4. Физические условия формирования земной атмо- 
сферы определили очевидный кусочный характер зависимости Т(п) . Общеизвестно, что в тропо- 


сфере температура линейно падает с высотой на ^^ —6,5 К/км вплоть до границы стратосферы 
(^^11 км). Далее до ^20 км она практически не меняется и держится на уровне -^216,65 К 
(—-56,5 °С) до ^20 км. Далее следует практически линейный рост до ^30 км и затем рост до 
^.50 км по криволинейной зависимости. 

Табличная форма задания параметров позволяет представлять данные с любой требуемой 
точностью, но не всегда удобна. При использовании таких данных в компьютерных расчётах 
(например, для планирования траекторий движения стратосферных летательных аппаратов и 
управления ими) удобнее иметь зависимости атмосферных параметров от высоты и других фак- 
торов в виде аналитических функций. Это избавляет вычислитель от необходимости работы с 
громоздкими массивами и решения задач интерполяции. 

В связи с этим может быть поставлена задача построения математической модели стан- 
дартной атмосферной температурной зависимости от высоты по точечным табличным данным, 
проиллюстрированным на рис. 4. При реализации первого этапа метода сиЕо/ие выделены четыре 
участка ломаной кривой: 0-11; 11-20; 20-30 и 30-45 км. Для каждого участка методами регрес- 
сионного анализа найдены аппроксимирующие их функции 


Т, (1) =288,13-6,4887 . 1; 

Т, (№) = 216,65; 

Т, (в) = 196,751 + 0,992 -1; 

Т, (№) = 492,914 21,748. +0,6165. 12 —5,1413.13. 


(12) 


Графики этих функций построены на рис. 4. Хорошо видно, что они практически совпадают с экс- 
периментальными точками на своих участках. Поэтому из них могут быть построены ИВФ с помо- 
щью МВФ, координатные параметры которых определены границами интервалов: №, =0; А =11; 


№, =20; В, =30; №, =45. Точностный параметр = =0,001 подобран экспериментально. 


* Гипотетическое вертикальное распределение температуры, давления и плотности воздуха в атмосфере Земли, которое 
по международному соглашению представляет среднегодовое и среднеширотное состояние [19]. 
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В результате из функций Т, (х) мультипликативно по формуле (4) строятся ИВФ, а из них 


по формуле (3) — итоговая «склеенная» функция: 
т (В) =Т, (1). ^(В, В, В, =)+Т, (В)-^(В, В, В, =) +Т. (В) - (В, В, В, =)-+Т. (В) -^(В, В, В, =). (13) 
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аппроксимации участков 
р ЕЙ 


0 5103 1х10* 1.5х10% 2х10* 2.5х10\ 3х10* 3.510 4х10' Й,м 


Рис. 4. Международный температурный стандарт атмосферы (МСА) и функции ф, (х) аппроксимации его четырёх участков 


График этой непрерывно-кусочной функции построен на рис. 5 и визуально полностью 
совпадает с точками табличных значений температуры. Численная проверка погрешности ап- 
проксимации всех табличных данных международного стандарта атмосферы показала, что по- 
грешность полученного математического описания не превышает 0,03 К. Этот показатель значи- 
тельно превышает точность оценки, необходимую для инженерно-технических и даже научных 
расчётов с использованием предоставляемого стандартом параметра. Можно построить матема- 
тические модели и других стандартизированных параметров атмосферы, что сделает возможным 
введение в стандарт аналитических их представлений наряду с табличными. При необходимости 
точность описания табличных данных можно повысить, уточнив границы участков. 
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Рис. 5. Результат сиЕд/ие аппроксимации температурного стандарта атмосферы 


Выводы. Изложен вполне регулярный метод сиЁд/ие аппроксимации, отличающийся существен- 
ной новизной решения задачи. Известные методы аппроксимации строятся лишь на приближении 
некоторой функции к описываемой точечной совокупности или её участку. Такой подход имеет 
слабые стороны. При использовании для всей совокупности единого выражения теряется точ- 
ность описания. При использовании множества моделей для каждого аппроксимируемого участка 
возрастает громоздкость и сложность описания. Сложность кусочного описания состоит в необхо- 
димости согласования моделей по краевым условиям. 

Новизна метода сиЁд/ие аппроксимации заключается в следующем. Точностные, согласу- 
ющие и объединяющие функции разделяются между: 
— математическими моделями описания выделенных интервалов ф, (х), мультипликативно вы- 


деляющими чисто интервальные модели # (х) выражениями ^(х, 5 


1 


— аддитивной формой объединения моделей # (х) в единую функцию / (х). 


р 
Это разделение, а также применение мультипликативного выделения фрагмента аппрок- 
симирующей функции существенно упрощает решение общей задачи и открывает практически 
неограниченные возможности для повышения точности аппроксимации. 
Метод не ограничен возможностями только одномерного применения. Совершенно оче- 
видно, что он вполне применим и в многомерном случае, когда составляющие функции многих 
переменных аддитивно разделены по аргументам. 
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В рассмотренном примере метод продемонстрировал компактность, простоту реализации и 
хорошие точностные показатели и вполне может быть рекомендован для использования при ап- 
проксимационном математическом описании экспериментальных зависимостей очевидного кусоч- 
ного типа, в особенности при наличии резких сломов между участками. 
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АРРКОХТМАТТМС МАТНЕМАТТСАЕ МОШЕ ОЕМЕГОРМЕМТ АССОКОТМС ТО РОТМТ 
ЕХРЕК1МЕМТАЕ ОАТА ТНВОУСН “СОТ-СТЦЕ” МЕТНОО” 


В. А. Меудоге 


А 5о/ивоп {ю те ргоМет оп аезсн тд ехрептетеа!у офатеа аерепсдепсе$з /5 соп!егед. Тве аштог5 тетпоа 15 
Базе ироп де@та зоте /оса/ арргохйтаНоп$ оЁ падтеп оР {Тезе геайоп5, апа те! абаШме гедисйоп Ю а 
по апа/уНса! ехргез$юп. Тиз ейесеЕ 15 аЕегттеа итд зреса! ‘а!осайпд” Кипсвоп$ ИтИтад пе дота! о! поп- 
его дей оп Гог еасй о! {те арргохтайоп ипсйопз. Тре теёоа 15 са!еа “сиЕ-дие” ассогАта №0 {те арр/еа 
рипар/ез. ТВе созезЕ апа/одие о! {те ргорожей тетод 15 5р/те арргохтаНоп. Ноиеуег те “сиЕ-дие” тетпоа #5 
тисй тоге адараЫЕ, ас /# 15 Бопаеа {0 пейпег {пе питьВЕГ о! <р/пе-арргохита Ме ротёз, пог {пе Кипсвоп огаег 
арргохитайпд Ё1е агеа$. ТПе огаег о! те роупоти! арргохтапЕ, ог апо{ Тег арргохтайпд ГипсНоп, аз ие! а$ #5 
5ргисиге Гог еаср хе, сап Бе агЬйгагу. АпоПег адуатадеоиз атегепсе оЁ "сиЕ-дие” арргохтайоп соп$Ё5 т а 
по апауйс поайоп оЁ те и/по! р/есеилсе Кпсйоп т$еаа оЁ аейптд а уесог р/те-ГипсНоп Егоидв а сит- 
Бегсоте эузет о! едиаНопз. Т№/5 ейесЕ ваз Бееп асШеуеа итд ап апа/у#са! Рипсйоп арргохитайпа апа рага- 
тесайу агЬйгат/ арргоасйтд {пе Неаизае $ер Кипсоп. ТВе апауйса! ап питепса! зНис!ез ог {те ргорегиез 
апа {те ебесЁ5 ог! аррутд йе ргоро5ей те{То4 аге гезийед. Тре оМатеа гезий5 аге /изёгаёеа ий" пе зресйс 
{есйгиса/ затр/е арр/сайоп$ о! {пе теёроа {ю ргасёса!/ ргоМетз, {абшаг апа дгар!иса! аа. 

Кеуигога&: ехрептета! аерепсЕепсе, р/есеил5е ипсйоп, арргохтайоп, тир/савийу, аа ит, оегепнарКу, 
апа/уйс Гипсйоп, рагатеёс арргоасн. 


* Тне могК {5 Чопе оп {Не {ПетаНс гезеагсН р!ап оЁ {не Митеу оЁ ЕдисаНоп ап@ 54епсе оЁ 1е Виз$ап ЕедегаНоп по. 2.3.13 
“Ва!апсеЯ затр!пд те{о4 Гог этшаНоп то4е!та ргоетз ог дупатгс ргосез5ез п аЬщед обес.” 
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